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ABSTRACT 

Observations have shown that the mass of nitrogen dioxide decreased at both southern and northern 
midlatitudes in the year following the eruption of Mt. Pinatubo, indicating that the volcanic aerosol had 
enhanced nitrogen dioxide depletion via heterogeneous chemistry. In contrast, the observed ozone response 
showed a northern midlatitude decrease and a small southern midlatitude increase. Previous simulations that 
included an enhancement of heterogeneous chemistry by the volcanic aerosol but no other effect of this 
aerosol produce ozone decreases in both hemispheres, contrary to observations. The authors’ simulations 
show that the heating due to the volcanic aerosol enhanced both the tropical upwelling and Southern 
Hemisphere extratropical downwelling. This enhanced extratropical downwelling, combined with the time of 
the eruption relative to the phase of the Brewer-Dobson circulation, increased Southern Hemisphere ozone 
via advection, counteracting the ozone depletion due to heterogeneous chemistry on the Pinatubo aerosol. 


1. Introduction 

The 15 June 1991 volcanic eruption of Mount Pina- 
tubo injected about 20 Tg of sulfur dioxide (S0 2 ) into 
the stratosphere (Bluth et al. 1992), up to an altitude 
of about 30 km (McCormick and Veiga 1992). The 
S0 2 transformed into about 30 Tg of sulfate aerosol 
(McCormick et al. 1995), increasing the stratospheric 
aerosol loading by orders of magnitude over back- 
ground. This perturbation decreased with time but a 
signature of the eruption was detectable in the atmo- 
sphere for about 5 yr (Robock et al. 2000). 

Aerosol from Mt. Pinatubo reached both the North- 
ern Hemisphere (NH) and Southern Hemisphere (SH), 
changing the stratospheric chemistry and dynamics. The 
volcanic sulfate provided additional surfaces for the 
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heterogeneous chemical reaction that transform N 2 0 5 
into HN0 3 , resulting in the depletion of nitrogen dioxide 
(N0 2 ) (e.g., Fahey et al. 1993; Tie and Brasseur 1995; 
Portmann et al. 1996; Solomon et al. 1996; Tabazadeh 
et al. 2002). This heterogeneous reaction generally leads 
to an increase in CIO and OH, and subsequently to an 
increase of ozone depletion (e.g., Tie and Brasseur 1995). 

Additionally, heating by the volcanic aerosol changed 
the radiative balance and dynamics of the atmosphere 
(e.g., Robock and Mao 1992; Graf et al. 1993; Stenchikov 
et al. 2002; Robock et al. 2007; Karpechko et al. 2010), 
intensifying the upwelling in the tropics and downwelling 
in the extratropics (Kinne et al. 1992; Pitari and Mancini 
2002; Aquila et al. 2012). 

The observed depletion of stratospheric N0 2 in both 
hemispheres during the years following the eruption 
(Johnston et al. 1992; Van Roozendael et al. 1997; 
Danilin et al. 1999) provided evidence that the volcanic 
aerosol had enhanced the heterogeneous chemistry at 
all latitudes. Surprisingly, the observed ozone response 
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Fig. 1 . Zonal mean of the ozone column anomaly (%) from 
TOMS observations. The baseline is the 1979-90 column ozone 
mean from TOMS. Effects due to the seasonal cycle, QBO, ENSO, 
solar cycle, and changing chlorine and bromine have been removed 
by using a seasonally varying regression analysis (Stolarski et al. 
1991). The green and purple triangles mark the locations and times 
of the eruption of Mt. Pinatubo and Cerro Hudson, respectively. 

to the volcanic perturbation was different between the 
NH and SH. While column ozone generally decreased in 
the NH, an increase of the ozone column was detected at 
southern mid- and high latitudes during the year follow- 
ing the eruption (Zerefos et al. 1994; Randel et al. 1995; 
Randel and Wu 1996). Figure 1 shows the anomaly of the 
total ozone column due to the eruption of Mt. Pinatubo 
as derived from Total Ozone Mapping Spectrometer 
(TOMS) observations. 

Several model studies attributed the observed NH ozone 
depletion to the enhancement of heterogeneous chemistry 
because of the volcanic aerosols (e.g., Tie and Brasseur 
1995) but could not explain why the same heterogeneous 
chemistry did not affect the SH ozone concentration 
(Stolarski et al. 2006). Some studies suggested explanations 
other than a chemistry perturbation for the asymmetry of 
the ozone response. Randel and Wu (1996) showed that the 
quasi-biennial oscillation (QBO) increased ozone in the 
extratropical SH during the 1991/92 winter, but the effect is 
not large enough to explain the observed increase. Fleming 
et al. (2007) and Telford et al. (2009) successfully simulated 
the ozone behavior using observed meteorological fields. 
These studies attribute the absence of ozone depletion to 
interannual dynamic variability, which masked the Pina- 
tubo aerosol chemical effect in the SH. However, these 
studies cannot distinguish between natural interannual 
variability and circulation changes forced by the volcanic 
perturbation. Poberaj et al. (2011) performed a multiple 
linear regression analysis to the Chemical and Dynamical 
Influences on Decadal Ozone Change (CANDIDOZ), 
stating that volcanically induced chemical ozone depletion 
was overcompensated by the QBO and by a pronounced 
Eliassen-Palm (EP) flux anomaly. 

Here, we use a free-running global chemistry-climate 
model (section 2) to separate the photochemical and 


dynamical contributions to the ozone and N0 2 anomalies 
induced by the volcanic perturbation alone (section 3). 

2. Model and simulation 

All simulations are performed with the Goddard 
Earth Observing System Chemistry-Climate Model 
(GEOSCCM). The general circulation model is the 
Goddard Earth Observing System, version 5 (GEOS-5; 
Rienecker et al. 2008), a system of component models 
integrated using the Earth System Modeling Framework 
(ESMF). For these simulations, GEOS-5 is coupled to 
the Georgia Institute of Technology-Goddard Global 
Ozone Chemistry Aerosol Radiation and Transport 
(GOCART) module (Colarco et al. 2010) and a strato- 
spheric chemistry module (Pawson et al. 2008). The 
resolution is 2.0° latitude X 2.5° longitude with 72 ver- 
tical layers from surface to 0.01 hPa (—80 km). This 
version of GEOSCCM does not simulate the QBO. 
The model is forced with observed sea surface tem- 
peratures and sea ice concentrations (Reynolds et al. 
2002). Aquila et al. (2012) includes a more detailed 
description of the model and an evaluation of the 
simulation of the Mt. Pinatubo cloud’s transport within 
GEOSCCM. 

GOCART computes the transformation of S0 2 into 
sulfate aerosol. The simulations shown here use a pre- 
scribed aerosol surface area density (SAD) from the 
Stratospheric Aerosol Measurement II (SAM II), 
Stratospheric Aerosol and Gas Experiment (SAGE), 
and SAGE-II data (Eyring et al. 2008) for heterogeneous 
chemistry. GEOSCCM also includes an option for cal- 
culating the aerosol SAD online using the simulated mass 
of sulfate aerosol and the relative humidity. 

We simulate the Mt. Pinatubo eruption by injecting 20 
Tg of S0 2 between 16 and 18 km in the volcano’s model 
grid box on 15 June 1991. No other aerosol sources are 
included in the simulation. We performed four model 
experiments, each composed of 10 simulations with 
different sets of initial conditions typical of the year 
2000, when the values of effective equivalent chlorine 
were similar to the early 1990s (Clerbaux and Cunnold 
2007). The results shown here are ensemble averages. 

The first experiment (experiment REF) is a control 
ensemble that does not include any volcanic perturba- 
tion. This experiment uses prescribed aerosol SAD from 
SAM II and SAGE observations from 1979 for the 
heterogeneous chemistry — a period when the strato- 
spheric aerosol layer is relatively unperturbed (Thomason 
et al. 1997). The heterogeneous chemistry in experiment 
REF is representative of the chemistry that takes place 
on the stratospheric background aerosol in absence of 
major volcanic eruptions. 
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Table 1. List of performed experiments. 


Experiment 

Radiatively interactive 
aerosol 

Year(s) for sulfate area 
density 

Perturbation to the 
chemistry 

Perturbation to the 
dynamics 

REF 

No 

1979 

No 

No 

CHEM 

No 

1991-95 

Yes 

No 

DYN 

Yes 

1979 

No 

Yes 

FULL 

Yes 

1991-95 

Yes 

Yes 


The second experiment (experiment DYN) includes 
radiatively interactive aerosol; that is, the aerosol can 
modify the simulated meteorology. As in experiment 
REF, experiment DYN uses aerosol SAD from 1979. 
Hence, this experiment perturbs the stratospheric dy- 
namics but does not perturb the stratospheric chemistry. 

The third experiment (experiment CHEM) does not 
include radiatively interactive aerosol and uses aerosol 
SAD from SAGE-II data appropriate for the simulated 
year; that is, the heterogeneous chemistry takes place on 
an aerosol layer perturbed by the Mt. Pinatubo eruption. 
This experiment includes a Pinatubo perturbation of the 
stratospheric chemistry but no perturbation of the 
stratospheric dynamics. 

The fourth experiment (experiment FULL) includes 
radiatively interactive aerosols and aerosol SAD for the 
Pinatubo period. This experiment includes the volcanic 
perturbation to both the stratospheric dynamics and 
chemistry. 

The comparison of the experiment FULL with the 
control experiment REF identifies the complete per- 
turbation of ozone and N0 2 due to Mt. Pinatubo. 
Comparisons of experiments DYN and CHEM with 
REF isolate the anomalies of ozone and N0 2 due to the 
effects of the eruption on the atmospheric dynamics and 
heterogeneous chemistry, respectively. Table 1 sum- 
marizes the performed model experiments. 

3. Results 

Figure 2 shows the observed deviations of the strato- 
spheric N0 2 column (black line) over Lauder, New 
Zealand, from a visible and ultraviolet (UV/Vis) spec- 
trophotometer (Johnston and McKenzie 1984). The 
baseline for the calculation of the anomalies is the ob- 
served monthly means over the years 1997-2003, which 
is a period free of major volcanic eruptions. The plot 
shows the simulated zonal mean anomalies at 45°S of the 
stratospheric N0 2 column, calculated as the difference 
between the experiments FULL (red), CHEM (blue), 
and DYN (green) and the experiment REF. The DYN 
experiment including only the perturbation to the dy- 
namics does not show any significant perturbation of 
N0 2 . Both the CHEM (blue) and FULL (red) experiments 


produce a decrease of stratospheric N0 2 in good agreement 
with the observations. The perturbation of N0 2 is domi- 
nated by the volcanic effect on the stratospheric chemistry, 
which is relatively symmetric between the hemispheres 
(Fig. 3). Heterogeneous reactions on sulfate aerosols 
explain the repartitioning between NO x (NO + N0 2 ) 
and the reactive nitrogen reservoir NO^ (NO + N0 2 + 

no 3 + 2N 2 0 5 + hno 3 + ho 2 no 2 + ciono 2 + 

BrON0 2 ) that is observed after the eruption of 
Mt. Pinatubo at northern and southern midlatitudes 
(Fahey et al. 1993 and Koike et al. 1994). 

The simulated N0 2 depletion is smaller in magnitude 
and peaks 1 month later than the observations. Koike 
et al. (1994) presented results from the Lawrence Liver- 
more National Laboratory (LLNL) (Wuebbles et al. 1991) 
and Atmospheric and Environmental Research (AER) 
(Rodriguez et al. 1991, 1994) models, which also used 
aerosol fields from SAGE-II and also underestimate the 
N0 2 depletion. We performed experiments using the on- 
line calculation of the aerosol SAD (not shown), and 
found a larger and earlier depletion of N0 2 in better 
agreement with the observations. 

GEOSCCM produces a reasonably good simulation 
of the vertical distribution of the Pinatubo aerosol in the 
posteruption period (Aquila et al. 2012). Figure 4 (left 
panel) displays the SAGE-II 525-nm zonal mean 



Fig. 2. Stratospheric N0 2 column anomalies (%) at 45°S. The 
black line shows the anomalies at Lauder, New Zealand (45.5°S), 
calculated from UV/Vis observations. The baseline is a 1997-2003 
monthly climatology calculated from these Lauder observations. 
The solid black line is significantly different from zero at ler level. 
The simulated 45°S anomalies for FULL - REF (red), DYN - 
REF (green), and CHEM - REF (blue) are also shown. Shading 
shows the standard deviation of each ensemble. The solid red, 
green, and blue lines are significantly different from zero at the 
95% confidence level, while the broken lines are not significant. 
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Fig. 3. Zonal mean stratospheric N0 2 column anomalies in ex- 
periment FULL with respect to experiment REF. The gray shaded 
areas are not significantly different from zero at 95% confidence 
level. 


extinction coefficients for the 1991 September-November 
(SON) period. Features to note include the large extinc- 
tions in the tropics and the SH midlatitude profile with 
double peaks at about 40 and 200 hPa. The online model 
simulation produces a somewhat different SAD vertical 
profile in the SH midlatitudes, with peaks at 20 and 
200 hPa (Fig. 4, right panel, blue line). The upper peak 
is slightly higher than is observed in SAGE-II. The 
analysis of SAGE-II data by Kinnison et al. (1994) also 
revealed this two-peak SH vertical structure. The pre- 
scribed monthly mean profile from SAGE-II and used 
herein (red line; Eyring et al. 2008) only exhibits a single 
peak in aerosol SAD at about 70 hPa. The reasons for 
the discrepancy between this SAGE-II single peak 
profile shape (Fig. 4, right panel, red line) and the 
SAGE-II extinction observations (left panel) possibly 
results from the coarse averaging resolution used in 
deriving the prescribed values indicated by the red line. 

The volcanic perturbation to the dynamics modifies 
the ozone concentrations. Figure 5 compares the 30°-60°N 
and 30°-60°S zonal total ozone column anomalies from 
TOMS data (Herman et al. 1996; McPeters et al. 1996) with 
the GEOSCCM simulations. The observed anomalies 
(black) are calculated as the deviation from the ozone 
column mean over 1979-90 after removing the effects of 
the seasonal cycle, QBO, El Nino-Southern Oscillation 
(ENSO), solar cycle, and changing stratospheric chlo- 
rine and bromine. The data show a negative NH 
anomaly (Fig. 5, top panel) from December 1991 that is 
significantly different from zero (all solid lines in the 
figure are significantly different from zero) for about 
2 yr after the eruption. The experiment CHEM also 
shows a negative NH anomaly, while the experiment 
DYN has no significant ozone anomaly. The negative 
anomaly present in the experiment FULL is not signif- 
icantly different from the anomaly of CHEM at the 95% 
confidence level. Hence, the simulated NH anomaly is 
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Fig. 4. (left) Zonal mean vertical profile of the SAGE-II aerosol 
extinction coefficient at 525 nm. (right) Zonal mean vertical pro- 
files at 30°-60°S of the prescribed aerosol surface area density 
adopted for our simulations (red line, Eyring et al. 2008) and the 
online simulation from GEOSCCM (blue line). The blue shaded 
areas show the standard deviation of the 10-member ensemble 
used for the online calculation of the SAD. Both panels are 
September, October, and November 1991 averages. 

primarily caused by the Pinatubo-forced heterogeneous 
chemical effect. 

In the SH, TOMS recorded a positive anomaly for about 
1 yr after the eruption (Fig. 5, bottom panel). This SH 
positive anomaly is simulated by both the DYN and FULL 
experiments. By June 1992 the DYN positive ozone 
anomaly is gone, while CHEM produces a significant 
negative perturbation starting about March 1992. This 
CHEM ozone response is similar in both hemispheres. The 
ozone anomaly in FULL is essentially the sum of the DYN 
and CHEM anomalies. The ozone anomaly in FULL is 
significantly different from CHEM at 95% level during 
two periods: October 1991-January 1992 and March-July 
1992 (gray shaded areas of Fig. 5). These differences 
demonstrate that the ozone anomaly from the Pinatubo 
dynamical response is larger or comparable to the chem- 
ical perturbation for about 1 yr posteruption. 

The positive ozone anomaly at SH midlatitudes in ex- 
periments DYN and FULL is induced by the dynamical 
perturbation. The absorption of largely longwave radia- 
tion by the volcanic aerosol leads to enhanced upwelling 
in the tropical lower stratosphere. By continuity, this 
enhanced tropical upwelling enhances the downwelling in 
the extratropics (Aquila et al. 2012). The ozone anomaly 
due to advection by this enhanced circulation is negative 
in the tropics and positive at midlatitudes (Fig. 6). This 
volcanic perturbation of the dynamics lasts as long as 
there is a strong gradient of the aerosol concentration 
between the tropics and midlatitudes. This mean circu- 
lation perturbation by Pinatubo was created by the 
volcanic aerosol radiative forcing (Aquila et al. 2012). 
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Fig. 5. Zonal mean total ozone column anomalies for (top) 30°- 
60°N and (bottom) 30°-60°S. The black line shows TOMS anom- 
alies (%) after removing the effects of the seasonal cycle, QBO, 
ENSO, solar cycle, and changing chlorine and bromine (see Fig. 1 
for details). The solid black lines are significantly different from 
zero at ler level. The simulated anomalies for FULL - REF (red), 
CHEM - REF (blue), and DYN - REF (green) are also shown. The 
color shaded areas show the standard deviation of each ensemble. 
Solid red, blue, and green lines are significantly different from zero at 
the 95% confidence level. The gray shaded times show where the 
FULL and CHEM simulations are significantly different from one 
another at the 95% confidence level. The green triangle marks the 
month of the Pinatubo eruption. 


The anomalous ozone distributions (Fig. 7) are shown 
in the immediate posteruption period [June- August 
(JJA) 1991, left panel] and the late posteruption period 
(SON 1992, right panel) by differencing the FULL and 
REF simulations. In JJA 1991 the perturbation is 
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Fig. 6. Zonal mean anomalies of total ozone column between the 
FULL and REF simulations. The gray shaded areas are not sig- 
nificantly different from zero at 95% confidence level. 

completely dominated by the dynamics response, which 
we depict with the streamlines of the residual circulation 
anomaly. The increased tropical upwelling lifts air with 
lower ozone concentration, creating a negative equato- 
rial anomaly centered at 20 hPa or 24 km. At the same 
time, the downwelling south of the equator creates an 
SH positive ozone anomaly. This positive anomaly is 
located in the SH because of the phase of the Brewer- 
Dobson circulation at the time of the eruption, which is 
directed toward the winter (SH) hemisphere. We per- 
formed an experiment initiating a Pinatubo-like erup- 
tion on 15 January 1991, as described in Aquila et al. 
(2012). There, the positive ozone anomaly appeared in 
the NH (not shown) — compatible with the different 
phase of the Brewer-Dobson circulation. 

In SON 1992 (Fig. 7, right panel) the simulated ozone 
anomaly is mainly due to the perturbation to the 
chemistry. In the middle stratosphere the ozone con- 
centration is increased, because the presence of the 
volcanic aerosol suppresses the NO x cycle that destroys 
ozone. In the lower stratosphere the presence of the 
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Fig. 7. Altitude vs latitude distributions of the zonal mean ozone anomalies (%) in (left) immediate post-Pinatubo 
period of JJA 1991 and (right) late post-Pinatubo period of SON 1992. The anomalies are calculated by differencing 
the FULL and REF simulations. The streamlines in the left panel show the residual circulation anomalies. 
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volcanic aerosol enhances the HO x and C10 x cycles and 
causes ozone depletion, as was also described in Tie and 
Brasseur (1995). 

4. Conclusions 

The lack of observed ozone depletion due to the 
eruption of Mt. Pinatubo in the Southern Hemisphere, 
in spite of the clear depletion of N0 2 , has been an 
outstanding puzzle for many years (Douglass and 
Fioletov 2011). We have shown that the perturbation 
of the stratospheric dynamics by the eruption of 
Mt. Pinatubo is responsible for the lack of an observed 
ozone decrease in the SH during the first year after 
the eruption. In the SH the dynamical response to the 
volcanic perturbation dominates the changes in ozone 
column during the first 6 months after the eruption and 
fades away starting about January 1992. The chemical 
response, instead, produces significant changes in the 
SH starting about 1 yr after the eruption. The pertur- 
bations to the chemistry and to the dynamics have an 
additive effect, resulting in the lack of ozone depletion 
at southern midlatitudes in the year following the 
eruption. 

The N0 2 anomaly is completely driven by the volcanic- 
induced chemical perturbation and is insensitive to the 
dynamic perturbation. The reason is the much shorter 
time scale of the heterogeneous chemistry for de- 
pleting N0 2 compared to the time scale for depleting 
ozone, together with the weak N0 2 vertical gradient, 
such that changes of vertical advection do not induce 
large perturbations. 
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